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Abstract 

We have experimentally investigated field induced aggregation of nonmagnetic particles confined 
in a magnetic fluid layer when rotating magnetic fields were applied. After application of a magnetic 
field rotating in the plane of the fluid layer, the single particles start to form two-dimensional (2D) 
clusters, like doublets, triangels, and more complex structures. These clusters aggregated again and 
again to form bigger clusters. During this nonequilibrium process, a broad range of cluster sizes was 
formed, and the scaling exponents, z and z' , of the number of clusters N(t) ~ t^'and average cluster 
size S(t) ~ t z were calculated. The process could be characterized as diffusion limited cluster- 
cluster aggregation. We have found that all sizes of clusters that occured during an experiment, 
fall on a single curve as the dynamic scaling theory predicts. Hovewer, the characteristic scaling 
exponents z' , z and crossover exponents A were not universal. A particle tracking method was 
used to find the dependence of the diffusion coefficients D s on cluster size s. The cluster motions 
show features of Brownian motion. The average diffusion coefficients < D s > depend on the 
cluster sizes s as a power law < D s >cx s 7 where values of 7 as different as 7 = —0.62 ± 0.19 and 
7 = —2.08 ± 0.51 were found in two of the experiments. 
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I. INTRODUCTION 



Colloidal aggregation phenomena are intersting subjects of study for both theoretical 
and technological reasons. In systems with short range interactions the main aggregation 
features are well understood (author?) [1]. Diffusion limited cluster-cluster aggregation 
(CCA) model (author?) [2, 3] and dynamic scaling theory (author?) [4] explain well the 
scaling properties during aggregation. It was found that these models, initialy developed 
for systems with short range interactions, can be used in systems where dipole-dipole in- 
teraction is dominant, for example aggregation of magnetic microspheres (author?) |s[|f3], 
aggregation of nanoparticles in magnetic fluid (author?) jj] and aggregation of magnetic 
holes (author?) |8|]. These experimental results show scaling of the significant parameters 
and features typical of CCA. On the other hand, the corresponding exponents may deviate 
slighly from the known models, and the reasons for this is still not understood. 

Our previous results (author?) jg] served as motivation for this study. Aggregation of 
magnetic holes in constant magnetic fields was interpreted in frames of the CCA model and 
dynamic scaling theory. The scaling exponent z ~ 0.42 for the cluster size dependence S(t) ~ 
t z was found for particles of diameters d — 1.9 and 4 /mi. This value of the exponent is slighly 
lower than exponent values predicted by theory (z = 0.5) (author?) [2j or exponents found 
by computer simulations (z = 0.5 for isotropic and z = 0.61 for anisotropic aggregation) 
(author?) [9]. Under certain experimental conditions (i.e., particles with larger diameter 
d = 14 fim) the exponent z was close to or lower than the exponent value corresponding 
to a transition from 2D aggregation to ID aggregation [z = 1/3). Based on these optical 
observations we know that for a constant magnetic field clusters move in 2D but they grow 
only in one dimension. 

The constant magnetic fields induced formation of long chains of particles (author?) 
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To determine correct scaling exponents we may take into account hydrodynamic 
effects (author?) We proposed to apply rotating magnetic fields to ensure quasi isotropic 
properties inside the magnetic fluid (MF) sample. In this case clusters can move in 2D space 
and can grow as 2D compact objects and thus the hydrodynamic correction is less important. 

The dynamic properties of a few magnetic holes (author?) [lOj in rotating magnetic fields 
show interesting phenomena, for example nonlinear response of bound pairs of magnetic holes 
(author?) [11], complex braid dynamics (author?) 12|, and equilibrium configurations 
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FIG. 1: (colour online) Experimental setup used to study aggregation of magnetic holes. A rotating 
magnetic field H(t) is applied in the plane of the magnetic fluid. 



of rotating particles without contact between particles (author?) [13J. In a precessing 
magnetic field paramagnetic particles dispersed in a drop of water self-assemble into two- 
dimensional viscoelastic small clusters (author?) 141 ]. 

In the present study field induced aggregation of many magnetic holes has been observed. 
In Sec. [Ill the experimental equipment and the methods used are described. The Sec. ITTTI 
deals with the results concerning the determination of the scaling exponents and charac- 
terization of the diffusion behaviour of individual clusters by tracking of their motions. In 
Sec. IIVI we summarize the general features and try to explain the non-universal scaling 
exponents. Our conclusions follow in Sec. [V] 



II. MICROSCOPIC OBSERVATIONS 

The experimental setup shown in Fig. [T] consists of an optical microscope (Nikon Op- 
tiphot), two pairs of coils, and a carefully prepared thin layer sample. Alternating curents 
were supplied to the coils in order to produce a magnetic field rotating in the horizontal 
plane of the sample. Microscopic observations were captured by a CCD camera (Q-Imaging 
Micropublisher 5) with resulution 2560 x 1920 pixels. 

The sample size was about 20 x 20 mm 2 . A layer of magnetic fluid of thickness ap- 
proximately 50 fim was confined between two glass plates and sealed. The kerosene based 
magnetic fluid (author?) ]15] had the following physical properies: density p = 1020 kgm 3 , 
susceptibility \ = 0-8, saturation magnetization M s = 20mT, and viscosity rj = 6x 10 3 Nsm 2 . 
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Monodisperse polystyrene microspheres of diameter d = 3/zm were dispersed in the MF layer 
in order to create magnetic holes in presence of magnetic fields. 

Without a magnetic field the particles are homogeneously dispersed in the layer and they 
can move freely. After some time a very low fraction of particles may randomly join to other 
particles and a few dimers were observed (author?) [10]. Their volume fraction is very low 
in comparison with the volume fraction of single particles. However, before the application 
of the rotating magnetic field, a short magnetic field pulse perpendicular to the sample 
(coils are not shown in Fig. [T]) was applied in order to destroy these dimers and to create 
a monodisperse initial size distribution of particles. This initial stage of the experiment is 
not shown in Fig. El 

The rotating magnetic field H(£) = (H x , H y ) within the x — y plane had the compo- 
nents: H x = Hsin(ut) and H y = H sin(u;t + vr/2). The amplitude of the magnetic field was 
constant H = 793 Am -1 , and angular velocity u = 251 s _1 . The temperature during the ex- 
periments was T = 293 K. The effective volume susceptibility including the demagnetization 
correction for spherical magnetic holes was Yeff — x/(l + 2x/3) = 0.63. The dimensionless 
interaction strength parameter (author?) [8( was A = 90. Here, A = U^ x /kT, where U^ x 
is maximal dipolar energy of two joined dipolar particles, k is Boltzmann's constant and 
T is the temperature. Thus, dipole-dipole interaction among magnetic holes was dominant 
over the thermal fluctuations. 

The rotating magnetic field induced an aggregation of the microspheres. The process took 
place via the joining of single particles into dimers, trimers and formation of 2D clusters 
consisting of many particles. These new clusters aggregate again and formed bigger clusters. 
A typical aggregation dynamics is shown in Fig. [2j A few samples with approximately the 
same layer thickness were investigated. The volume fractions of particles were low, in the 
range = 0.0014 - 0.0064. 

In order to analyze the digital images, a C programing language code and open graphical 
libraries were used. Several thousands of digital pictures have been analyzed in a distributed 
manner in a computational grid. We have analyzed the motion of individual particles during 
aggregation using our own tracking algorithms written in the Python programing language. 
The main advantage of the algorithm that was used is that it can track positions of new 
clusters which are results of the aggregation. 
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FIG. 2: Optical micrographs of aggregation of nonmagnetic microspheres with diameter d = 3/mi 
in magnetic fluid at different times: (a) t = 600 s, (b) t = 3000 s, and (c) t = 3400 s after the 
magnetic field was switched on. The applied rotating magnetic field had an amplitude |H| = 793 
Am -1 and frequency / = 40 Hz. The optical view covers a sample area of about 368 /mix 274 /mi. 
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III. RESULTS 



Microspheres inside a magnetic fluid layer witout magnetic field behave as nonmagnetic 
particles dispersed in fluid. They perform random Brownian motion. In this case, aggre- 
gation events are rare due to the low particle concentration. Thus, in the initial stage of 
the experiments the microsperes are homogeneously dispersed in the layer of MF and the 
cluster size distribution is unimodal. 

After application of the external magnetic field the microspheres begin to behave as 
interacting magnetic holes. They have induced magnetic moments which are oppositely 
oriented to the external magnetic field. When the energy of dipole-dipole interaction among 
two arbitrary spheres is larger than the thermal energy of the spheres, as quantified by the 
dimensionless interaction strength A = 90 in the present case, field induced aggregation 
starts. 

During the aggregation complex motions of microsperes and clusters consisting of many 
microspheres were observed. Clusters containing regularly ordered particles were formed 
and small irregular clusters relatively quickly relax to highly ordered structures. Based on 
the optical observation the complex modes of motion of microsperes and clusters may be 
classified as; i) joining of two clusters together followed by a very slow relaxation of the 
microsperes in the new cluster into a more ordered structure; ii) extremely slowly swivelling 
of all clusters in the same direction as the rotating magnetic field, followed by packing into a 
compact disk form; and iii) small random motions of the clusters induced by random forces 
resulting from interactions with the local cluster environment. 

We have observed that clusters of all sizes can join together and form bigger cluster 
which is the basic feature of cluster-cluster aggregation. The cluster-cluster aggregation 
model (author?) [2] predicts the scaling properties of the total number of clusters N(t) and 
mean cluster size S(t). The total number of clusters is defined as N(t) = ^2 s n s (t) where 
n s (t) is number of clusters of size s at time t. The mean cluster size S(t) is defined as: 

S(t) _ Es n s(t)s 2 

where s is cluster size. In our case the cluster size s is given by the number of particles 
which belong to the cluster. 

The aggregation process in Fig. [2] was studied in more details. In Fig. [3] (a) we can 
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FIG. 3: (color online) (a) The total number of clusters N(t) and the mean (weight average) cluster 
size S(t) [in units of number of spheres] versus time, (b) The scaling function g(x) = s 2 n s (t) 
obtained from the cluster size distributions n s {t) during the time interval t = 200 — 3400 s. 



see that the total number of clusters N(t) and mean cluster size S(t) show power law 
dependencies N(t) ~ t z ' and S(t) ~ t z . The power-law behaviour was found for the time 
interval t = 200 — 3400 s. The scaling exponents were determined as z' = 0.65 ± 0.01 and 
z = 0.64 ±0.01. 

Based on dynamic scaling theory (author?) [4] all number of clusters n s (t) observed 
during aggregation can be scaled into a single, universal curve or scaling function g(x) 
denned as: 



n s {t) ~ s- 2 g(s/t z ). 



(2) 



It is expected that g(x) ~ x A for i < 1. 



All the cluster number curves n s (t) during the time interval t = 200 — 3400 s fall onto the 
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FIG. 4: (color online) (a) The total number of clusters N(t) and the mean (weight average) cluster 
size S(t) [in units of number of spheres] versus time, (b) The universal scaling function g(x) ~ x A 
(x = s/S(t), x < 1) calculated for the time intervall t = 200 — 80000 s. 

single curve shown in Fig. [3j From that the characteristic scaling exponent A = 1.44 ±0.08 
was found. 

We measured several samples, however, other samples behaved in a different manner. The 
scaling exponents z, z', and crossover exponent A were different from the results presented 
above. Typical results for a sample that shows a different type of behaviour are shown in 
Fig. H](a). Here the scaling exponents were found to be z' = 0.40 ±0.03 and z = 0.34 ±0.02. 
Similarily to the case discussed above, the cluster numbers n s (t) (t = 200—80000 s) that were 
measured for this sample could be scaled onto a single curve as shown in Fig. H] (b), but the 
scaling exponent A was nearly twice as large is in the former case, A = 2.75±0.06. Also in 
this case the visible dynamic behaviour was consistent with diffusion limited cluster-cluster 
aggregation but with clearly different scaling exponents from those above. 

The results presented in Figs. [3] and 0] show that the scaling exponents for this system 
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FIG. 5: (color online) Tracks of particles and clusters in the experiment shown in Figs. [2] and [3] 
during the time interval t = — 3400 s. Red lines belong to short tracks that are shorter than 10 
time steps (the time step At = 30 s). Blue tracks are longer than 10 time steps. The sample area 
is about 368 ^mx274 /im. 

can not be universal. In order to understand this unexpected result we have investigated 
the motions of individual clusters in more detail. 

The complex motion of a cluster was simplified by considering only the motion of its 
central mass point. There are effects that can change the position of the central mass point 
with nearly no motion of the cluster as a whole. For example, after joining of two clusters a 
rearrangement of particles in the new cluster (see the case i) discussed above) takes place. 
We assume that these disturbing changes are smaller than the influence of random local 
forces (case hi)) that essentially contribute to the cluster motions. A very slow rotation of 
a cluster (case ii)) does not change the position of the central mass point. 

Cluster tracks shown in Fig. (Fig. [U]) were determined for experimental data shown 
in Fig. [3] (Fig. HI). We see in Figs. [5] and [6] that the clusters moved in two directions, the 
tracks are complex and show features of Brownian motion as expected. 

For Brownian particles it is characteristic that their motion is well described by 

< |r| 2 >oc (3) 

where r is the distance vector between an initial position and the position after time t and 
D is the diffusion coefficient. We checked the validity of Eq. [3] for the cluster tracks and 
determined the relation D s oc< |r| 2 > ft for any cluster of the size s. At each experiment 
we analyzed about 1000 tracks and found that this equation is valid with a cluster-size 




FIG. 6: (color online) Tracks of particles and clusters during aggregation in the time interval 
t = — 3400 s. Red lines belong to short tracks that are shorter than 10 time steps (the time step 
At = 30 s). Blue tracks are longer than 10 time steps. The sample area is about 368 ^mx274 //m. 

dependent diffusion coefficient D s . 

We found that D s clearly depended on the cluster size s as shown in Fig. [7J For the case 
in Fig. [7J b) the values of the diffusion coefficient D s fall in a broad range covering nearly 
three decades. These ranges of values of D s are significantly larger than possible errors of 
measurement. An average diffusion coefficient < D s > was calculated for any cluster size s 
and the results (Fig. [7J were fitted to a scaling law < D s >oc s 7 . For the data presented in 
Figs. [S] and E] the diffusion scaling exponent 7 = — 0.62±0.19 was found. On the other hand, 
for the data in Figs. HJandElthe diffusion scaling exponent is clearly higher 7 = — 2.08±0.51. 



IV. DISCUSSION 



We have observed that application of a rotating magnetic field on a 2D magnetic many 
holes system causes field induced aggregation. The clusters move and can grow in both 
dimensions, which is different from the case of a constant magnetic field where clusters 
are free to move in both dimensions but only grow in one dimension as determined by the 
external magnetic field. The results show that the system is in a nonequilibrium state and its 



t), develope and 
2| as was shown 



characteristic quantities, number of clusters N(t) and average cluster size S 
show scaling according to a cluster-cluster aggregation model (author?) 
in Figs. EK) andHH). Both the scaling properties and the broad cluster size distributions 
found, as well as the existence of a scaling function g(x), are main signatures of cluster- 
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FIG. 7: (color online) The diffusion coefficients vs. cluster size s. (a) For the experimental 
data shown in Figs. and the diffusion coefficient follows < D s >oc s -°- 62±019 . (b) For the 
experimental data shown in Fig. 0] and [6] the diffusion coefficient scales as < D s >cx s ~ 2 - 08±0 - 51 . 

cluster aggregation and dynamic scaling theory (author?) ^]. 

In many cases the cluster-cluster aggregation mechanism leads to formation of complex, 
fractal-like objects (author?) [1]. However, in the present case the structure of the aggre- 
gates are simpler with a compact internal organization. In relatively strong magnetic fields, 
small clusters are regular 2D objects with well ordered structure of the particles inside the 
clusters. Extremely slow cluster rotations and rearrangement of particles inside the new, 
bigger clusters have been observed. As a consequence of these effects the clusters are packed 
into regular objects with a nearly close-packed, triangular structure of the spheres. The 
cluster diffusion coefficients do not depend on the direction in which the clusters move, i.e., 
hydrodynamic corrections are not important as they are in case of constant magnetic fields 
(author?) fl. 
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In the basic cluster-cluster aggregation model (author?) [2J] it is assumed that the dif- 
fusion coefficient is 7 = —1 and the corresponding scaling exponent is z = 0.5. We have 
determined two distinct values for the diffusion coefficients 7 and scaling exponents z as 
a result of the two clearly different types of behaviour observed in our experiments. The 
relationship between 7 and z for both values of 7 follows the equation z = 1/(1 — 7) which 
has been found in other aggregation models. For 7 = —0.62 ± 0.19 (7 = —2.08 ± 0.51) 
we computed z = 0.62 [z = 0.32). These scaling exponents agree well with exponents z 
determined directly from the time dependence of S(t), z = 0.64 ± 0.01 and z = 0.34±0.02, 
respectively. 

Unfortunately, at present we are not able to explain why similar experiments on approxi- 
mately the same samples (concentrations, layer thickness etc.) show scaling exponents with 
values that come in two clearly separated ranges and diffusion exponents 7 which are differ- 
ent from the expected value 7 = — 1. Thus, the scaling exponent z is either clearly lower or 
higher than the theoretically predicted value z = 0.5. 

In an earlier study of a similar system of magnetic holes in a constant magnetic fields 
(author?) [8] it was found that for small microspheres (diameters d = 1.9 — 4.0/zm and 
interaction strength A = 8 — 370) the scaling exponents z and z' were approximately equal 
z = z'and typically slightly lower than 0.5: 0.38< z, z' < 0.54. However, for larger particles, 
d = 14/im (A = 1040 — 10600) the values of z and z' increased with the value of the 
dimensionless interaction strength A from about 0.1 to 0.6, and correspondingly the value 
of A decreased from above 3.0 to ~ 1.5. Thus, depending on the particle size the scaling 
exponents changed from being nearly constant and near the theoretically expected values 
to being strongly non-universal. Although the present particles are within the diameter 
and A ranges which showed nearly universal behaviour in Ref. (author?) jg], the magnetic 
interactions are very different (anisotropic in the former and isotropic in the present) and 
then the range of A for which the behaviour in non-universal, seems to be changed. It is 
unclear why the diffusion conditions, as quantified by the values of the diffusion coefficients 
7, were so different in the two typical cases reported here. It may possibly be related to fine 
details in the interaction between the microspheres and the glass plates confining the system. 
In principle the magnetic holes should be repelled from the confining walls (author?) [h]] 
but if for some unknown reason a small fraction of the particles become attracted or even 
loosely attached to the walls, this would slow down the diffusion as shown by the anomalous 
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value 7 — 2 in one of the analyzed cases. Extremely small values of < D s > could indicate 
that some of the particles are trapped in the sample volume or on the sample glass boundary. 

V. CONCLUSIONS 

Diffusion limited cluster-cluster aggregation of magnetic holes has been induced by a 
rotating magnetic field. The main features of the experimental results are well described 
by a diffusion limited cluster-cluster aggregation model and dynamic scaling theory. The 
experimental conditions were designed in effort to have a well defined model of a low con- 
centrated many body system where long range interactions are dominant. At present, the 
reason why two main aggregation regimes were observed is not clear. This resulted in scal- 
ing exponent values clearly different from those predicted by theory for systems with short 
range interactions. This difference in behaviour was further confirmed by unusual values 
of the scaling exponent 7 of the diffusion coefficient that were found by cluster tracking. 
An open question remains: Why do the isotropic, long-range particle-particle interactions 
suppress the diffusion regime where the size-dependence of the diffusion coefficient scales 
with 7 = — 1? This will hopefully be clarified in future studies. 
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